Emergent Z 2 gauge symmetry and spin-charge separation in one dimensional physics. 
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We demonstrate that in the scaling limit the phenomenon of spin-charge separation as encoun- 
tered in Luttinger liquids defined on the lattice can be associated with the emergence of local Ising 
symmetry. This Z2 gauge field is of geometrical nature (sublattice parity) and the Luttinger liquid 
corresponds with the critical point where the symmetry is on the verge of becoming local. We argue 
that in a true one dimensional superconductor the local symmetry will be fully realized as long as 
the geometrical structure associated with spin-charge separation is maintained. 
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It is quite well understood how to break symmetry, but 
a different matter is how to make symmetry. With the ex- 
ception of electromagnetic U(l), the manifest symmetries 
governing condensed matter physics are of a global na- 
ture. Nevertheless, the idea that gauge theory is behind 
the peculiarities of high T c superconductivity is quite in- 
fluential, especially so in the form of the slave theories 
of spin-charge separation |l],|2| . Here we will demonstrate 
that behind the phenomenon of spin-charge separation 
in one spatial dimension a dynamically generated local 
symmetry of the Ising (Z2) kind is hidden. This gauge 
principle is different from those envisaged in the slave 
theories ||. It is rooted in geometry: the gauge fields 
parametrize the fluctuations of the geometric property 
bipartiteness. A lattice is called bipartite when it can be 
subdivided in two sublattices A and B, such that all sites 
on the A sublattice are neighbored by B sublattice sites 
and vice versa. This division can be done in two ways 

( A-B-A-B • • • and B-A-B-A ■ ■ •) defining 

a Z2 valued quantity 'sublattice parity', p = ±1. Bcthe- 
Ansatz results show that in the scaling limit sublattice 
parity acquires the status of separate dynamical degree 
of freedom. The quanta of electrical charge are bound 
to flips in the sublattice parity which in turn alter the 
space seen by the spin system. Sublattice parity turns 
into a Z2 gauge field when the quantization of charge 
is destroyed by true superconducting long range order, 
while the Luttinger liquid is right at the quantum phase 
transition where this local Z2 symmetry emerges. 

The above claims are based on a re-analysis of exact 
Bethe-Ansatz results for the one dimensional Hubbard 
model |||. Some time ago, Ogata and Shiba || discov- 
ered a simple but most peculiar property of the Bethe 
Ansatz wavefunction in the special case that U is very 
large (U — > 00 limit). They showed that the wave func- 
tion factorizes in a charge part ipSFi depending on where 
the electrons are, and a spin part tpu which is merely 
depending on the way the spins are distributed, 
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The charge part ipSF is nothing else than the wave- 
function of a non-interacting spinless-fermion system 
where the coordinates xt refer to the actual positions of 
the electrons/spinless fermions. The spin part ijjjj is iden- 
tical to the wavefunction of a chain of Heisenberg spins 
interacting via a nearest-neighbor antifcrromagnetic ex- 
change. Although this factorization property gives a 
precise meaning to the notion of spin-charge separation, 
charge and spin are actually not quite independent from 
each other. In ipu only the positions of the up spins 
are needed and these correspond with the coordinates jji. 
The surprise is now that the coordinates yt do not refer 
the original Hubbard chain, but instead to a new space: 
a lattice with sites at coordinates x±,X2, .-.jXjv given by 
the positions of the charges in a configuration with am- 
plitude 4>SF- 

The above gives away that the quantum dynamics of 
interacting electrons generates a geometrical structure 
analogous to the fabric of general relativity. Let us visu- 
alize this for a representative example (Fig. 1). Consider 
N electrons on a chain with L sites under the condition 
that N < L such that the charge configurations can be 
specified by the locations of the holes. A charge config- 
uration in the full Hubbard chain ('external space'), has 
an amplitude ipsF in the wavefunction with the coordi- 
nates of the dots corresponding with the a^'s. The spin 
system sees a different 'internal space' obtained from the 
full space by removing the holes together with the sites 
where the holes are located, substituting the hole and its 
site with an antifcrromagnetic exchange between the sites 
neighboring the hole (the 'squeezed space' ||). 

As in general relativity, physics is derived from re- 
lations between the different reference frames, but the 
geometry involved in the Ogata-Shiba case is obviously 
much simpler than the geometry of fundamental space- 
time. This simplification makes it possible to parametrize 
matters in terms of a simple gauge theory. In which 
regards are the full chain and the squeezed chain differ- 
ent? The squeezed chain is obviously shorter than the full 
chain and this is a simple dilation: a distance x measured 
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in the full chain becomes a distance px in the squeezed 
chain (p — N/L, the electron density) when the 
lattice constant. The other aspect is also simple, but 
less trivial in its consequences. The spin system is a 
quantum-antifcrromagnet and it is as such sensitive to 
the geometrical property of bipartiteness. Consider what 
happens with the sublattice parity, as introduced in the 
first paragraph. For the Heisenberg spin chain, as the 
one in squeezed space, a redefinition of p = 1 — 1 does 
not carry any consequence ('pure gauge'). However, sub- 
lattice parity becomes alive in the mapping of squeezed 
space into full space (Fig. 1). 'Fix the gauge' in squeezed 
space by choosing a particular sublattice parity, and con- 
sider what happens when it is unsqueezed. The holes are 
inserted, and because every hole is attached to one site, 
every time a hole is passed the sublattice parity flips. 

The above is true for every instantaneous charge con- 
figuration. However, the ground state is a superposition 
of many of these configurations. In 1+1D charge can- 
not break translational invariance and since this fluctu- 
ating charge is 'attached' to the sublattice parity flips, 
the full space which is observable by external observers 
(experimentalists) should be considered as a fluctuating 
geometry. However, this is a very simple fluctuating ge- 
ometry because all what is fluctuating is the property 
of bipartiteness. This geometrical fluctuation can be 
parametrized in terms of a field theory controlled by the 
simplest of all local symmetries: Z2 gauge theory 
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FIG. 1. The emergence of the Z2 'sublattice parity field' 
in the geometrical 'squeezing' operation. The open circles re- 
fer to some configuration of holes. The black dots define the 
embedding space of the spin system in both the full- (upper) 
and squeezed (lower) lattice. The unsqueezing operation can 
be parametrized in terms of binding of the electric charge 
quantum to flips in the Z2 valued sublattice parity field. 

Let us start out on the geometrical level and ask the 
question if a correlation function can be defined acting 
on the full Hubbard chain which can measure the 'true' 
spin correlations associated with squeezed space. Since 
all what matters is sublattice parity this can be achieved 
by simply multiplying the spin operator by a factor — 1 
every time a hole is passed, thereby removing the sublat- 
tice parity flips from the spin correlations. Define stag- 
gered magnetization as M(xi) = (— l) Xi S(xi) where S 
is the spin operator (S z = n-f — ri|,S l+ = ejej) and the 
charge operator n Xi = n It | + n Xi i taking the values 0, 
1 and 2 for an empty-, singly- and doubly occupied site, 



respectively. The correlation function we are looking for 
is §, 

O^fa - Xj \) = (M\ Xl )^-l +l e^-^M*{ Xj )). (2) 

The operator exp(i7r(l — n(xi)) takes the value +1 for a 
singly occupied spin/squeezed space site while it is —1 
for a charge (hole, or doubly occupied) site. By multi- 
plying these values on the interval X4 < xi < Xj all the 
minus signs associated with the sublattice parity flips are 
removed from the spin correlations. 

Although the 'string' operator ~ nexp(i7r(l — ri)) is 
non-local it can be evaluated straightforwardly using the 
techniques introduced by Parola and Sorella Q. It is 
easily shown that 

(^(oh-i^P 1 "" ^ 1 s'(xi)) 

Xi+l 

= £ PfM*;)(-i) Xi+1 ~^W^ - 1) (3) 

where Oueis. is the spin correlator of the Heisenberg 
chain, while Pg F (xj) = ( n(0) n{ Xi ) £Q^| =0 n(xi) - 
Xj)}sf is the probability of finding j spinless fermions 
in the interval [0, x*]. This factor causes the additional 
decay of the spin correlations due to the charge fluctua- 
tions in the standard spin correlator. However, it is easily 
shown that it is precisely compensated in O top by the fac- 
tor (— i^i+i-zj coming from the string operator and we 
find the result, asymptotically exact for large distances 
x, 

O top (x) = ^/ 2 {px). (4) 

Does this reveal the spin-only problem in squeezed space? 
We have to correct for the fact that the density of 
staggered spin in full space is reduced by a factor p 
as compared to squeezed space while in addition dis- 
tances are measured by x/p instead of x. After rescaling, 

Otop led { x ) = jhr 1 ^ 2 ^) which is indeed the behavior of 
the staggered spin correlation function of a Heisenberg 
spin chain at large distances. 

Let us now consider the well known || asymptotic be- 
havior of the 'normal' staggered spin correlations for the 
Hubbard chain ||, 

O a {\x i -x j \) = {\M*{x i )M'{x j )\) (5) 
~ C °1t ( £ )x) . (6) 

X + " 

What matters is that the algebraic decay is more rapid 
than for the correlator Eq. (^. The charge-stiffness K p 
is associated with the decay of the charge correlations, 
(|n(xj)n(0)|) ~ cos(2ex)/x K p, and K p > for all den- 
sities p / 1. Hence, the charge fluctuations influence 
the spin correlations because the charge is attached to 
the sublattice parity flips, and they do so by a simple 
multiplicative factor 1 / x p . 
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Is the above an accident of the strongly coupled 
case? Let us consider the case where the interac- 
tions are completely vanishing: the spinfull fermion 
gas on the lattice. It is straightforward to show 

V** -1 [1— n(x ■)] 

that in this limit ( S z ( Xl ) (-l)^*^ 1 S Z (Q)} = 

( ~ 1) 4 X '' 1 [D{xi - 2)D{xi) ~ D(x 2 - l) 2 ] where D(x - 

i) = { (— 1)^-^3=0 *■ This determinant can be eval- 
uated numerically and we find that O t0 p(x) ~ 1/x. As 
®m ~ l/x Kp+K " ~ 1/a; 2 , while if p = 1 in the non- 
interacting case, it follows that there is no difference in 
the asymptotics between the U — > oo and the U = 
case and by continuity it has to be that the squeeze is 
true at least for all U > 0. The exact limiting results 
(C7 = 0, — » oo) are further fortified by bosonization for 
arbitrary p which once again yields Ot op ~ x~ Ks [|To|. It 
is remarkable that free fermions code for a phenomenon 
as involved as the Ogata-Shiba squeeze! 

We have now gathered the information needed to make 
the case that a local Ising symmetry is at work. The key 
is that the non-local charge string ~ IIexp(i7r(l — n)) can 
be rigorously identified with the Wilson line of a Z 2 gauge 
theory ||, with the logical implication that the different 
decay rates of Otop and are manifestations of the 
local symmetry. The charge sector can be parametrized 
for large U in terms of spinless fermions, and these can 
in turn be expressed in terms of Pauli matrices a a using 
a Jordan- Wigner transformation, 
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where /J. creates a spinless fermion, while their number 
n Xi = fl fxi is identical to the electron charge operator 
n defined in the above. Hence, by operator identity we 
can write the vertex operator as, 



r(X-n Xi ) = e »(7r/2)(l+<) = ff 3 



(8) 



As charge and spin are decoupled, we may shift the 
charge degrees of freedom to the lattice formed by the 
midpoints of the bonds between the sites of the origi- 
nal lattice (the bond-, or dual lattice, with coordinates 
Xi,Xi+i) while concurrently retaining the spin degrees of 
freedom on the original lattice. After these operations, 
our correlator Eq. (ji|) becomes, 

O top {\ Xl -x \) = (|M^[n <Xi;ai , >er < i!C( jM|.|). (9) 

This is nothing else than the gauge invariant correlation 
function of a theory containing a SU (2) 'matter' field M 
coupled to an Ising gauge field. The factor Hr 3 corre- 
sponds with the Wilson-line of the gauge theory (T is 
a line connecting Xi and Xj), which has to be inserted 
to keep the matter- field correlator gauge invariant. By 



the same token, the 'normal' spin correlator O^ is not 
gauge invariant, because the Wilson line is missing, and 
our observation that is more rapidly decaying than 
Otop signals the presence of local Z 2 symmetry. 

Let us explain this in some more detail Q. Imagine 
that an effective, long wavelength theory is realized with 
a Hamiltonian invariant under the following Z 2 gauge 
transformation (cr XitXi+1 is the value of the Z 2 variable 
on bond Xi, 



&Xi,Xi+-L * ^Xi^Xi^Xi-^-i^Jxi-^-i 

M Xi -» rj Xi M Xi 



(10) 



with arbitrary r\ Xi = ±1. Given this symmetry, some 
very general statements follow regarding the behavior of 
correlation functions. First consider the spin correlator 

Oa(\xi-Xj\) = (M*(xi)M*(xi)). S 
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under the transformation Eq. ( |10[ ) at either Xi or Xj it 
has to vanish: it is not gauge invariant. Since the gauge 
symmetry is emergent an energy scale should exist below 
which it becomes active. An energy scale implies a length 
scale Igauge and therefore 0^(x) ~ exp(-x/; gQU g e ). The 
unique way to construct gauge invariant correlation func- 
tions is by inserting a closing Wilson line, Ilrcr 3 , as done 



in Otop (Eq. 



The minus sign introduced by the 



gauge transformation of e.g. M{xi) is compensated by 
the minus sign introduced by the simultaneous transfor- 
mation of the bond variable on the bond leaving site Xi 
in the direction of site Xj. However, a single cr 3 is not 
gauge invariant and to maintain overall gauge invariance 
the string ncr 3 , ending at Xj, needs to be inserted. The 
behavior of Otop will depend on the details of the dy- 
namics of the gauge theory. However, one statement is 
always valid: Otop will decay more slowly then O ^ when 
the gauge invariance is present, and this is what we found 
to be the case in the Luttinger liquid. 

Nonetheless, there is a complication. As we argued, 
an unavoidable consequence of a full realization of the 
Z 2 gauge symmetry is that correlation function violating 
the gauge invariance should decay exponentially faster 
than the gauge invariant ones, and in the Luttinger liquid 
we found that this difference is a mere algebraic factor. 
Algebraic behavior of correlation functions is associated 
with yet another symmetry: scale invariance. Since Otop 
shows a slower algebraic decay than O^ it has to be that 
the Luttinger liquid is right at the critical point where 
the Z 2 symmetry turns on, which is associated with the 
divergence of the gauge length: l gauge — * oo. 

This is less mysterious than it might sound. It is clari- 
fying to view this from the geometrical (squeezed space) 
perspective. The Z 2 local invariance just means that it is 
undetermined how many sublattice parity flips occur per 
unit length. Since these flips are attached to the charge 
quanta Z 2 invariance is established when charge is a fluc- 
tuating quantity. Z 2 turns into a global invariance in a 
state where charge is locally conserved: the Wigner (or 
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'holon') crystal - it is easily checked that Ot op and Om 
behave identically at long distances when the holes are 
assumed to be localized. In one dimension the tendency 
to crystallize is universal ||. However, this order always 
turns into algebraic order due to the admixing of the 
zero-modes and algebraic order is a form of criticality. 
This turns into criticality of the gauge sector due to the 
charge-sublattice parity binding. 

Hence, the system can be driven away from the phase 
transition to the phase obeying global Zi by breaking 
explicitly translational invariance. Is it possible to drive 
it in the other direction, into a phase where the local Zi 
invariance is truly realized? The answer is simple: add 
a field stabilizing true superconducting long range order 
in the charge sector. The charge fluctuations which are 
fundamental to a system with phase order are of pre- 
cisely the right kind to protect the local Z 2 symmetry. 
By adiabatic continuity, a system showing true supercon- 
ducting order can be thought of as composed of lattice 
bosons b'(xi) carrying a charge q. In addition, it is a 
requirement that this charged boson is attached to the 
sublattice parity flip . For long wavelength purposes, 
the charge part of the ground state wave-function of the 
phase-ordered state may be written as 



\tp q ) = Hcju + vb\xi)]\vac) 



(where, in the low density limit, |u| 3> \v\) such that the 
superconducting order parameter '5 = u*v. It follows 
that (Zr is the length of the Wilson line) , 



(iP g \U r e 



i7tb' (xi)b(xi) 



2\lr 



= (M - l«r) 



(12) 



where I 



l/ln(|u| 2 — \v\ 2 ). This in turn implies, 



gauge 



0$(\xi - Xj\) = e 's»"s= O top {\xi - Xj 



(13) 



proving that a full local Z2 invariance is generated. The 
key is that the spontaneous breaking of charge U(l) sym- 
metry is the same as the spontaneous making of local Z2 
symmetry, gauging the spin sector! 

Besides the clarification of the symmetries underlying 
one dimensional physics, does this new symmetry princi- 
ple buy us new physics in 1+1D in the sense of 'new states 
of matter'? This is disappointing, for no other reason 
than that gauge theories are not particularly interesting 
in 1+1D ||. The reason is that in one spatial dimen- 
sion the gauge field can only carry kinetic energy: using 
Hamiltonian language the only gauge invariant operator 
available to describe the dynamics of the gauge field is 
a li,x i+x = ( CT i,x i+1 +°x" i ,x j+1 )/ 2 - Hence, the ground state 
is always confining (eigenstate of cr 1 ), giving rise to ex- 
ponential decay of all spin-related correlation functions. 
This can be directly checked for our superconductor. In 
analogy with Eq. ([?]), a 1 ~ + b which acquires a vac- 
uum amplitude because (ipq\b* + b\ip q ) is non-zero in a 



superconductor. The situation is more interesting in two 
(and higher) spatial dimensions, because the gauge field 
can now carry also a potential energy. A variety of distin- 
guishable phases becomes possible [p 12,l3|. Elsewhere 
we will analyze the 2+1 dimensional generalization of the 
above, clarifying the meaning of the cuprate stripes and 
suggesting a novel way of viewing the high T c enigma. 
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